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Pressure- and temperature-induced structural phase transitions of CaFe 2 As 2 and 
BaFe 2 As 2 studied in the Hund's rule correlation picture 
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With the proposed Hund's rule correlation picture, i.e. the fluctuating Fe local moments with the 
As-bridged antiferromagnetic superexchange interactions, the exceptional collapsed tetragonal phase 
and related phase transitions observed in CaFe2As2 are well understood. With the same framework, 
a pressure-temperature phase diagram is predicted for BaFe2As2 as well, in which a paramagnetic 
tetragonal and a collinear antiferromagnetic orthorhombic structures to nonmagnetic tetragonal 
structure transitions take place around 4-8 GPa and 10-15 GPa respectively, and a nonmagnetic 
tetragonal to a nonmagnetic collapsed tetragonal structure transition takes place over 26 GPa. Our 
study helps better understand the complex correlation among crystal structure, magnetism, and 
electronic structure in pnictides, a precondition to understand the superconductivity in pnictides. 

PACS numbers: 74.70.Xa, 74.20.Pq, 74.25.Dw, 74.25.Bt 



Since the discovery of superconductivity in LaFeAsO 
by partial substitution of O with F atoms below 26K[T], 
intense studies have been devoted to physical properties 
of iron-based pnictides both experimentally and theoret- 
ically. There are, so far, four types of iron-based com- 
pounds reported, showing superconductivity after dop- 
ing or under high pressures, i.e. 1111-type i?eFeAsO (Re 
= rare earth) [T], 122-type AFe 2 As 2 (A=B&, Sr, or Ca) 
[21, Ill-type BFeAs (B = alkali metal) [3j, and 11-type 
tetragonal a-FeSe(Te) [4]. All these compounds share 
the same structural feature that there exist the robust 
tetrahedral layers where Fe atoms are tetragonally co- 
ordinated by pnicogen or chalcogen atoms and the su- 
perconduction pairing may happen. A universal finding 
is that all these compounds are in a collinear antiferro- 
magnetic (AFM) order below a tetragonal-orthorhombic 
structural transition temperature [5J [3] except for a-FeTe 
that is in a bi-collinear AFM order below a tetragonal- 
triclinic structural transition temperature [7HS]. 

The above finding brings out a serious issue regard- 
ing the mechanism behind the structural and antiferro- 
magnetic transitions, and their underlying relationship. 
There are basically two contradictive views upon this is- 
sue. The first one [10] is based on itinerant electron 
picture, which thinks that there are no local moments 
and the collinear antiferromagnetic order is induced by 
the Fermi surface nesting that is also responsible for the 
structural transition due to breaking the four-fold rota- 
tional symmetry. On the contrary, the second one is 
based on local moment picture. The J\-Ji Heisenberg 
model was, phenomcnologically |llj and from strong elec- 
tron correlation limit [12] respectively, proposed to ac- 
count for the issue. At the same time, we proposed [T3] 
the fluctuating Fe local moments with the As-bridged an- 
tiferromagnetic superexchange interactions as the driving 
force upon the two transitions, effectively described by 
the 3\-Ji Heisenberg model as well. 



Our proposal embodies the twofold meanings shown 
by the calculations [13 : (1) there are localized magnetic 
moments around Fe ions and embedded in itinerant elec- 
trons in real space; (2) it is those bands far from rather 
than nearby the Fermi energy that determine the mag- 
netic behavior of pnictides, namely the hybridization of 
Fe with the neighbor As atoms plays a substantial role. 
Here the formation of a local moment on Fe ion is mainly 
due to the strong Hund's rule coupling on Fe 3d-orbitals 
[13j . In this sense, our proposal can be considered as 
the Hund's rule correlation picture. We emphasize again 
that the Arsenic atoms play a substantial role in our 
physical picture. Subsequently we successfully predicted 
from the calculations [7J, based on this Hund's rule corre- 
lation picture, that the ground state of a-FeTe is in a bi- 
collinear antiferromagnetic order, which was confirmed 
by the later neutron scattering experiment [S]. 

At ambient pressure, like other pnictides, CaFe2As2 
undergoes a transition from a nonmagnetically ordered 
tetragonal (T) phase to a collinear AFM orthorhombic 
(O) phase below 170K [H] (T T ^o)- Nevertheless, the 
neutron diffraction measurements [TBI, I16j further found 
that, at 50K, CaFe2As2 experiences another exceptional 
transition from a collinear AFM orthorhombic phase to 
a new nonmagnetically ordered but collapsed tetragonal 
(cT) phase with a dramatic decrease in both the unit 
cell volume and the c/a ratio when the applied pressure 
is larger than 0.35 GPa (To_, c r). The corresponding 
phase diagram is schematically shown in Fig. 1(a). Thus 
how to understand this phase diagram and accompanied 
magnetism, especially the two nonmagnetically ordered 
tetragonal phases, becomes an outstanding problem. 

In this Letter, we show that the exceptional collapsed 
tetragonal phase and related structure transitions, ac- 
companied by the magnetic transitions, found experi- 
mentally in CaFe2As2 can be well understood within 
the proposed Hund's rule correlation picture. The im- 



portance of the fluctuating local moments in the struc- 
tural transitions is elucidated. Moreover, we predict a 
pressure- and temperature-induced phase diagram for 
BaFe 2 As2 in the same framework. 

The previous calculations [1?\ show that the ground 
state of CaFe2As2 is in an orthorhombic structure with 
a collinear AFM order at low temperature and low pres- 
sure, in an excellent agreement with the experiments. It 
was naturally proposed, from the Hund's rule correla- 
tion picture, that the high-temperature tetragonal phase 
of CaFe2As2 is a paramagnetic (PM) phase with local 
magnetic moments around Fe atoms randomly orientated 
due to thermal fluctuations, which is supported by the 
universal linear-temperature dependence of static mag- 
netic susceptibility in iron pnictides [TS] and our calcula- 
tion illustrated below. Here, we further propose that the 
pressure-induced exceptional collapsed tetragonal phase 
of CaFe2As2 is a nonmagnetic (NM) phase, in which the 
Hund's rule in Fe atoms is overcome so that the Fe mag- 
netic moments are quenched. The corresponding calcu- 
lations reported below confirm our proposal and well de- 
scribe the exceptional phase diagram of CaFe2As2- 

Our density functional theory (DFT) calculations were 
done using the general gradient approximation (GGA) 
for exchange-correlation potentials [19], the projector 
augmented wave method [5D] , and a plane wave basis set 
up to 600 cV as implemented in the Vienna ab-initio sim- 
ulation package|21|. The Gaussian broadening technique 
was used and a fc-mesh of 16 x 16 x 4 or 12 x 12 x 4 was 
adopted to sample the Brillouin zone of 1 x 1 or V2 x V2 
supercell, respectively. The shape of the lattice and the 
internal coordinates of all ions were fully optimized at a 
given cell volume. A series of volumes, differing in 0.5% 
in each dimension from one to another, were calculated 
for each state, as reported in Figs. [2]and[3J 

In calculation, it is very difficult to directly simulate 
a paramagnetic phase. Nevertheless, the checkerboard- 
AFM and paramagnetic phases share the following im- 
portant features: (1) local moments around Fe atoms, 
(2) zero net magnetic moments in a unit cell, and (3) 
the same magnetic symmetry. Moreover, the effect of 
spin fluctuations on static properties in the paramagnetic 
phase may be accounted for by the checkerboard-AFM 
state. Thus as long as we study the static properties 
rather than the low-energy dynamics, the checkerboard- 
AFM state is feasible to properly model the paramag- 
netic phase. Accordingly, the orthorhombic collinear 
AFM (c-AFM), tetragonal checkerboard-AFM namely 
Neel order (denoted as Neel-AFM), and tetragonal NM 
states are considered to represent the low-temperature 
orthorhombic (O) phase, the high-temperature tetrago- 
nal (T) phase, and the pressure-induced collapsed tetrag- 
onal (cT) phase in CaFe2As2, respectively. 

Phase transition from O to cT is so unusual that phase 
O collapses to phase cT under a small external pressure. 
Figure [2] shows the evolution of the total energy with the 
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FIG. 1: (a) Sketch of the phase diagram of CaFe2As2. Four 
transitions are marked: transition from tetragonal (T) phase 
to orthorhombic (O) phase when the temperature decreases 
below 170K at ambient pressure (Tt->o), in which a long- 
range collinear antiferromagnetic order is established; this 
magnetic order disintegrates in collapsed tetragonal (cT) 
phase by applying a pressure(To-» c r); phase T recovers from 
cT when increasing the temperature (T c t^t); and phase T 
transforms again to cT by applying a large pressure (Tt-i c t)- 

(b) Supercells of phases T or cT in 14/mmm symmetry, and 

(c) supercell of phase O in Fmmm symmetry. Green, pink, 
and red spheres represent Ca, As and Fe atoms, respectively. 
Black arrows denote the magnetic moments of Fe. The mag- 
netic moments are quenched in phase cT. 



unit cell volume (E — V curve) in the c-AFM and NM 
states. The measured structural parameters for phases 
O and cT at 50 K were well reproduced, within errors of 
1-2%, by the DFT calculations in the c-AFM and NM 
states, respectively, while the Neel-AFM state cannot re- 
produce. One sees that the volume decrease is accompa- 
nied by a loss of magnetic moment in the c-AFM state, 
as shown by the red numbers in Fig. [21 Nevertheless, in 
the transition region, the magnetic moment is still large, 
roughly 1.6-1.7 fiB, which indicates that the full sup- 
pression of magnetic moments, after the c-AFM to NM 
transition, cannot be simply ascribed to pressure-induced 
full derealization of localized electrons solely. 

Gibbs free energy G determines the stability of phases, 
which is defined as G = E - TS + pV, where E, T, 
5, p, and V stand for internal energy, temperature, en- 
tropy, pressure, and volume, respectively. The transi- 
tion pressure of To^cT can be thus worked out by p = 
{{Enm — E c _afm) — T(Snm — S c -afm))/{V c -afm — 
Vnm), whereas the entropy terms are negligible at ex- 
tremely low temperature, e.g. 4K. It turns out that the 
transition pressure from c-AFM to NM state is given by 
the slope of the tangent line between the NM and c-AFM 
curves in Fig. [2j i.e. 0.58 GPa at OK, which excellently 
agrees with the experimental result of 0.6 GPa at 4K[To]. 
In contrast, this value is smaller by one order of magni- 
tude than another calculated value of 5.25 GPa working 
with the Fermi surface nesting picture 22J. Meanwhile, a 
dramatic cell volume deduction of about 19 A 3 is found in 
our calculations, consistent with the experimental value 
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FIG. 2: E - V curves of the NM and c-AFM states modeling 
phases cT and O, respectively. Magnetic moment of each data 
point in phase O is indicated by a red number. A downshift 
of 34 meV of the NM curve is made according to —TAS c t,o 
at 50K. This predicts transition pressures of 0.58 GPa and 
0.32 GPa at OK and 50K, respectively. 
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FIG. 3: E - V curves of the NM (black solid), Neel-AFM 
(red dashed), and constrained Neel-AFM (blue dotted) states. 
Magnetic moment of each data point in Neel-AFM state is 
indicated by a red number. A downshift of 78 meV of the 
constrained Neel-AFM curve, corresponding to the entropy 
difference between phases T and cT at 300K, is shown in the 
inset, which gives a transition pressure of 2.3 GPa. 



of 14.23 A 3 at 50K. 

At a finite temperature, the entropy terms need to be 
considered. The entropy change from phase cT to O 
(AS c t,o) is estimated as 0.67 meV/K[23]. This lowers 
the G of phase cT by 34 meV at 50 K, as illustrated by 
the gray curve in Fig. [2] Following the same procedure, 
another tangent line in the figure gives a transition pres- 
sure of 0.32 GPa and a similar volume deduction, again 
excellently consistent with the measured pressure in be- 
tween 0.24 and 0.35 GPa at 50 K[15]. 

Phase transitions between T and cT - Figure [3] shows 
the E-V curves of the Neel-AFM and NM states. We 
find that the NM state has a lower energy than the Neel- 
AFM state at the equilibrium. The Neel-AFM state is, 
however, more stable than the NM state at larger vol- 



umes, for example, at a volume of 176.83 A 3 close to 
the experimental value of phase T at 0.63 GPa and 200 
K|15j. The corresponding local moment increases to 1.03 
Hb- And the nearest neighbors Fe-Fe separation in the 
Neel-AFM state is only 0.03 A shorter than the experi- 
mental value while 0.14 A in the NM state. 

There are two phase transitions defined with phases T 
and cT, i.e. transitions T c t^t an d Tt^cT as marked 
in Fig. [T] To well understand these two phases, another 
spin-polarized calculation with magnetic moments con- 
strained at 1.0 fiB was performed, as reported by the 
blue dotted curve in Fig. [3J Because of containing an 
extra freedom degree of magnetic moments, the param- 
agnetic tetragonal phase T represented by the Neel-AFM 
state has a larger entropy S than the collapsed tetrago- 
nal phase cT represented by the NM state. This lowers 
Gibbs free energy in phase T more than in phase cT when 
increasing temperature. This change is presented by the 
inset of Fig. [31 in which the blue curve (constrained Neel- 
AFM) goes lower at a finite temperature by —TAS. As 
a consequence, the system takes transition from phase 
cT to phase T when the entropy change overcomes the 
volume expansion by increasing temperature, meanwhile 
the magnetic moments recover. The entropy change is es- 
timated as ~ 0.26 meV/K between phases T and cT|23j. 
by which a transition pressure of 2.3 GPa for T^^r is 
derived from the S'-corrected E — V curve (E — TS versus 
V in FigJ3]inset) at T=300 K, reasonably consistent with 
the experimental value of 1.5 GPa|16j. 

Collapse mechanism of CaFeiAsi - As analyzed above, 
the Fe magnetic moment is induced by the Hund's rule 
coupling on the Fe 3<i-orbitals, which is about 0.6-0.8 
eV/Fe. In transition To-± c t, the Hund's rule coupling 
must be overcome and compensated from Gibbs free en- 
ergy gain since the magnetic moments are quenched. In- 
spection of the calculations shows that the spin degener- 
ation pushes down the energy levels raised by the Hund's 
rule coupling, and then strengthens the bonding between 
the Fe 3d orbitals and As 4p orbitals, which makes sub- 
stantial energy gain to the NM state. Besides, there is an 
energy of ~0.1 eV/Fe (denoted as Er) contributed from 
other sources to compensate the Hund's rule coupling, as 
suggested by the total energy difference between the non- 
collapsed NM state and the c-AFM state of CaFe 2 As2 in 
the equilibrium. There are such possible sources as: (1) 
entropy change effect of —TAS; (2) volume shrink ef- 
fect of pAV; and (3) formation of new bonds that make 
further energy gain. 

The entropy difference between the c-AFM and NM 
states is very small and can be negligible in very low 
temperatures. As reported above, the cell volume con- 
taining four Fe ions is dramatically reduced by 19A 3 after 
transition To-> c t from the c-AFM to the collapsed NM 
state. In order to compensate Er in Gibbs free energy to 
favor the collapsed phase through the volume shrink ef- 
fect of pAV alone, an applied pressure of 3.4 GPa would 
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FIG. 4: Sketch of the predicted phase diagram of BaFe2As2. 



be required. However, it is noticed that CaFe2As2 is so 
compact that the nearest As- As distance d,As-As hi the 
equilibrium, as shown in Fig. [lib), is already very close 
to the distance that can form an As- As covalent bond, i.e. 
~2.9 A. Actually, the As- As bond indeed forms in the col- 
lapsed phase [53] . The formation of such a new bond con- 
tributes an additional energy gain of ~0.1 eV/Fe, as es- 
timated from the total energy difference between the col- 
lapsed and non-collapsed NM states of CaFe2As2- Such 
an energy gain is fairly comparable with E^. It turns out 
that the c-AFM to the collapsed NM transition can take 
place at a low pressure that pushes the nearest As atoms 
close enough to form bonds, due to the small size of Ca 
atoms and its resulting short inter-layer As- As distance. 
The collapsed phase is thus closely related to the forma- 
tion of the As- As bonds. In most pnictides, the nearest 
As-As distance in equilibrium is much larger than the 
one in CaFe2As2, differences in phase diagrams are thus 
expected, as found in BaFe2As2 reported below. 

Prediction on phase diagram of BaFe^Asi - Within 
the Hund's rule correlation picture, we predict a phase 
diagram of BaFe2As2, as shown in Fig. [4j which looks 
similar in shape to that of CaFe2As2- However, the cal- 
culation shows that the realization of an effective As-As 
bond needs a huge pressure (over 20 GPa) due to a large 
dAs-As hi BaFe2As2. Such a large pressure already pro- 
vides a sufficient Gibbs free energy gain through pAV to 
compensate the Hund's rule coupling, resulting in a non- 
magnetic state, not necessarily with an extra energy gain 
from the As-As bonding in collapsed structures. Con- 
sequently, there is no c-AFM directly to collapsed NM 
transition. Instead, there are two transitions with in- 
creasing pressure, first from the c-AFM to a normal NM 
state, then further to a collapsed NM state. The corre- 
sponding transition pressures are predicted following the 
same way adopted for CaFc2As2. They are, as shown 
in Fig. El (1) the paramagnetic tetragonal to the NM 
tetragonal phase transition at 4.0-7.2 GPa depending on 
temperature T; (2) the c-AFM orthorhombic to the NM 
tetragonal phase transition around 12.8 GPa at very low 
T; (3) the NM tetragonal to the collapsed NM tetragonal 
phase transition at 26-30 GPa. 



In summary, with the Hund's rule correlation picture, 
the exceptional phase transitions accompanied with the 
magnetic structures and properties in CaFe2As2 are well 
understood, as demonstrated using the density functional 
theory calculations. The calculated transition pressures 
of 0.58 (0.32) GPa at (50) K for the transition from the 
collincar antifcrromagnetic orthorhombic phase to the 
nonmagnetic collapsed tetragonal phase, and 2.3 GPa at 
300 K for the transition from the paramagnetic tetrago- 
nal phase to the nonmagnetic collapsed tetragonal phase, 
are highly consistent with the experimental values 0.6 
(0.25-0.35) GPa at 4 (50) K and 1.5 GPa at 300 K, re- 
spectively. Different from CaFe2As2, a nonmagnetic non- 
collapsed tetragonal phase is predicted to appear in the 
phase diagram of BaFe2As2 due to the larger atomic size 
of Ba than that of Ca, separating the nonmagnetic col- 
lapsed tetragonal phase from the collinear antifcrromag- 
netic orthorhombic phase and paramagnetic tetragonal 
phase. The corresponding transition pressures are pre- 
dicted as well and reported in phase diagram Fig. [4] 
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